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We study the relaxations, surface energies, and work functions of low index metallic surfaces 
using pseudopotential plane-wave density-functional calculations within the generalized gradient 
approximation. We study here the (100), (110), and (111) surfaces of Al, Pd, Pt, and Au and the 
(0001) surface of Ti, chosen for their use as contact or lead materials in nanoscale devices. We 
consider clean, mostly non-reconstructed surfaces in the slab-supercell approximation. Particular 
attention is paid to the convergence of these quantities with respect to slab thickness; furthermore, 
different methodologies for the calculation of work functions and surfaces energies are compared. 
We find that the use of bulk references for calculations of surface energies and work functions can 
be detrimental to convergence unless numerical grids are closely matched, especially when surface 
relaxations are being considered. Our results and comparison show that calculated values often 
do not quantitatively match experimental values. This may be understandable for the surface 
relaxations and surface energies, where experimental values can have large error, but even for the 
work functions, neither local nor semi-local functionals emerge as an accurate choice for every case. 



PACS numbers: 68.47.De, 71.17.Mb, 71.15.Nc 

I. INTRODUCTION 

Electronic structure calculations are becoming more 
widely applied to increasingly complex and realistic ma- 
terials systems and devices, reaching well into the domain 
of nanotechnology with applications ranging from metal- 
molecule junctions, carbon-nanotube field effect transis- 
tors, and nanostructured metals or semiconductors. 1,2 
For such complex applications, characterizing the prop- 
erties of the elementary building blocks becomes of fun- 
damental importance. Consider nanoscale constructions 
such as self assembled monolayers, where the underly- 
ing metal work function (and Fermi energy) are crucial 
to determining interfacial phenomenal Thus, it becomes 
critical to have a clear understanding of the fundamental 
properties of the metal lead surfaces - starting from sur- 
face energies, structural relaxations, and work functions. 

It is also important to highlight the limitations of den- 
sity functional theory when considering even only the 
building blocks, in this case with regard to surface prop- 
erties of metals. In general local and semi-local approxi- 
mations for the exchange-correlation functional perform 
well for metals, however this is not the case for metal 
surfaces^ For this purpose we proceed here with a de- 
tailed study of the properties of clean metal surfaces of 
interest, as lead materials, to nanotechnological applica- 
tions, and discuss the issues involved in obtaining con- 
verged estimates using electronic-structure modeling - in 
this case using density-functional theory (DFT) and the 
slab-supercell approximation - and the accuracy of these 
estimates compared to experimental values. 

First, surface relaxations can arise from the creation 
of a new surface - i.e. cleaving of a bulk in two. This 
leads to a smoothing of the charge density at the new 
surface which causes a net force on the outermost surface 



layer of ions pointing into the bulk. 5 DFT calculations 
in the slab-supercell approximation, assuming a pristine 
unit cell of the unreconstructed surface, can be effectively 
used to study these surface layer relaxations. 

Second, the surface energy is the energy required to 
cleave an infinite crystal in two - i.e. the amount of 
energy required to create a new surface. This is a dif- 
ficult quantity to determine experimentally because it 
usually requires measuring surface tension at the melt- 
ing temperature of the metal£ Theoretical determination 
of this quantity is relatively easy and particularly use- 
ful in studies of the relative stability of different surface 
facets! 7 - However surface energy calculations within DFT 
are sensitive to numerical errors arising from differences 
in Brillouin-zone sampling. 8 Methodologies for avoiding 
these errors have been proposed in the literature 8 -"^ and 
are discussed in later sections of this study. 

Third, the work function is the minimum energy 
needed to remove an electron from the bulk of a mate- 
rial through a surface to a point outside the material. In 
practice, this is the energy required at K to remove an 
electron from the Fermi level of the metal to the vacuum 
potential;- 1 Calculations of work function using DFT em- 
ploy this definition and determine the Fermi energy and 
vacuum potential from calculations of the metals in slab- 
supercell geometries. However, work functions calculated 
with slab approximations are known to have a depen- 
dency on the thickness of the slab, thus further analysis 
is required to extract bulk metal work functions from slab 
approximation. This dependency is well documented in 
some cases and is attributed to finite size effects arising 
from classical electrostatic interactions or from quantum 
size effects (QSE) i 12 i 13 i 14 Methodology to lessen such ef- 
fects is available in the literature^ 

In this paper we calculate the relaxations, surface en- 
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ergies, and work functions of the (111), (100), and (110) 
surfaces of Al, Pd, Pt, and Au and the (0001) surface of 
Ti within DFT. We examine the convergence of all three 
quantities with respect to the thickness of the slabs, fol- 
lowing the surfaces of Pd as an example. Furthermore, 
we employ and compare a number of different methods 
of calculation with respect to surface energy and work 
function. Well-converged values of these quantities for 
all of the metals and surfaces of this study are tabulated 
and compared, when available, to experimental quan- 
tities and/or DFT-LDA quantities from the literature. 
Some of the LDA generated data from the literature, and 
compared to here, were calculated with surface Green's 
functions mcthodologies i 9 ' 16 ! 17 These methods have been 
shown to be highly accurate, as well as inherently free 
from finite-size effects due to choices of boundary condi- 
tions. However, we feel that the flexibility and popularity 
of plane-wave DFT methods for more complex systems 
highlights the necessity for comparison of methods and 
functionals with well converged slab-supercell calculated 
quantities. 



II. METHODOLOGY 

First-principles calculations within density functional 
theory are carried out using the PWscf code of 
the Quantum-ESPRESSO distribution^ The Perdew- 
Burke-Ernzerhof (PBE) functional within the general- 
ized gradient approximation (GGA) is usedJ^ For the 
Brillouin-zone integration, we use a Monkhorst-Pack set 
of special k-points, 20 and Marzari-Vanderbilt smearing 21 
with a broadening of 0.02 Ryd. 

Ultrasoft pseudopotentials (USPP) are used for Pd, Pt, 
and Au. For the case of Pd and Pt the USPPs used have 
been generated with the RRKJ3 scheme^ 2 , with 10 va- 
lence electrons each in 4d 9 5s 1 and 5d 9 6s 1 configurations, 
respectively. The Au USPP have also been generated 
with the RRKJ3 scheme with 11 valence electrons in a 
5d 10 6s 1 configuration^ For the case of Ti the USPP was 
generated using the Vanderbilt scheme 23 with 12 valence 
electrons in a 3s 2 3p 6 4s 2 3d 2 configuration. A normcon- 
serving pseudopotential is used for Al with electrons in 
a 3s 2 3p 1 configuration. The kinetic energy cutoff for the 
plane wave basis are 32 Ryd for the wave function and 
512 Ryd for the charge density for Al, Pd, Pt, and Au. 
For Ti these values are 40 Ryd and 640 Ryd, respectively. 

Surfaces were constructed using a supercell with a thin 
slab of metal separated from its periodic images by a 
layer of vacuum. The size of this region is such that 
there are always ~ 16A of vacuum between the surfaces. 
For all the slabs considered, each layer in the unit cell 
contains one in-equivalent atom. For (111) surfaces a 
hexagonal cell with a base defined by ao[110] and do [101] 
and an inter-layer spacing of and a stacking of AB- 
CABC is used, where ao is the equilibrium lattice pa- 
rameter. The other surfaces are similarly constructed; 
for (100) a tetragonal cell is used with a base formed 



by -^[110] and -^[110] where the inter-layer spacing 
is and ABAB stacking. The (110) surface is con- 
structed from an orthorhombic cell with a base defined by 
^|[110] and do [100] with an inter-layer spacing of 
and ABAB stacking. For the case of Ti the hexagonal 
cell is used with a ABAB stacking in the z-direction. 

The remaining methodological items are discussed in 
the following results section. 



III. RESULTS 
A. Bulk Properties 

We calculate the lattice parameters and bulk moduli 
for Al, Pd, Pt, Au, and Ti. The calculations are per- 
formed for the total energy of the bulk system for a range 
of lattice parameters a. The total energy data are fit 
with the Murnaghan equation of stated to obtain the 
bulk moduli. For the case of bulk Ti, a hexagonal cell is 
used with a two atom basis to construct the hexagonal 
close packed (HCP) structure. Total energy calculations 
are performed for a range of c/a ratios and each of them 
fit with the Murnaghan equation of state. The results 
are found in Table HI Here we see that the lattice param- 
eters of the FCC Al, Pd, Pt, and Au are over estimated 
(< 2%) by the use of the PBE-GGA exchange correlation 
functional. Similarly, for the bulk moduli the calculated 
values are underestimated with respect to the experimen- 
tal values. However, for the calculations of these prop- 
erties of Ti, PBE slightly overestimates the volume of 
the unit cell as well as over-estimating the bulk modulus. 
This type of behavior can be seen in a study of similar 
metals when calculated using the PBE functional and a 
full potential linear augmented plane wave (FLAPW) all 
electron methodology^ From the same FLAPW study, 
we see good agreement with our pseudopotential calcu- 
lations for the lattice parameters (within 0.03A) and the 
bulk moduli (within 5 GPA) of Al, Pd, and Pt. 10 The re- 
mainder of this work the equilibrium lattice parameters 
calculated with PBE, found in Table [l] are used. 

Recently in the literature, it has been shown that Pd 
calculated within GGA is incorrectly described as having 
a magnetic ground state in the bulk i 25 i 26 To this end we 
also have calculated the bulk and surface properties of 
Pd with spin-polarization. We also find that the ground 
state is described as magnetic, however, when applied 
to the surface properties, the introduction of magnetism 
has a small effect. These calculations are reported and 
discussed in the Appendix of this paper. 



B. Surface Relaxations 

Surface relaxation is characterized as the percent 
change AaV,- of the spacing between layers i and j versus 
the equilibrium layer spacing d,Q. For the (111), (100), 
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TABLE I: Calculated lattice parameters and bulk moduli of 
the metals considered in this study compared with experi- 
mental values. 
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FIG. 1: Layer relaxations for the top six layers of Pd(100) as 
a function of slab thickness. The numbers reported are the 
percent change, Adij, of the spacing between layers i and j 
versus the initial layer spacing do (for the (100) case do = 
a /V2). 



(110), and (0001) surfaces d is a /\/3, a Q /y/2, a y/2/^ 
and co/2, respectively. 

As an example, the layer relaxations for Pd(100) as a 
function of slab thickness can be seen in Fig. [TJ Con- 
vergence of the relaxed layer spacing is achieved with 
increasing number of layers, and a slab thickness of 13- 
layers assures all relaxations are converged to within 
0.1%. We show in Tabic HTl the top three surface layer 
relaxations for all surfaces for 13-layer slabs. As ex- 
pected, the relaxation of the surface layers is related to 
the density of packing,— with larger relaxation for the 
less-densely packed surfaces, with patterns of multilayer 
relaxation that become noticeable as we go from (111) 
to (100) to (110) surfaces. For the case of Ti(0001), 
the surface is hexagonal close-packed, however the large 
layer spacing leads to large relaxations. Comparison to 
low-energy electron diffraction (LEED) and other exper- 
imental values shows that our well converged values do 
not correspond directly to the experimentally observed 
values. However, regarding first layer relaxation or ex- 



pansion, the qualitative trend is captured in all but the 
(100) surfaces. 

Finally, we note that for this study we are typically 
concerned with clean non-reconstructed surfaces (larger 
supercells in the in-plane dimension could allow for sur- 
face reconstruction). Some of the more widely know sur- 
face reconstructions are the (22 x \/3)-reconstruction of 
the Au(lll) surface^ the (2 x 1) missmg-row recon- 
struction of the (110) surfaces of Pt and Au ; 33 ' 34 and 
a number of reconstructions of the (100) surfaces of Au 
and Pt, the most studied being the (100)-hex^& To 
this end we investigate for comparison the Au(110) and 
Pt(110)(2 x l)-missing row reconstruction for 13-layer 
slabs. We find that in both cases the first-layer relaxation 
more closely matches that of the experimental values in 
Table HU -19.7% and -18.6% for Au and Pt, respec- 
tively. However, no improvement is seen in the second 
layer relaxation. Finally, the surfaces of the 4d-metal 
Pd as well as the Ti(0001) surface are know to be less 
likely to reconstruct at low temperatures. For the case of 
surface relaxations, neither the LDA nor GGA exchange- 
correlation functionals appear to represent a more accu- 
rate choice when compared to experimental values. 

C. Surface Energies 

The surface energy is the energy required to create a 
new surface, and as mentioned earlier, it is a difficult 
quantity to determine experimentally In our calcula- 
tions the surface energy a can be determined by taking 
the energy difference between the total energy of a slab 
and an equivalent bulk reference amount, as seen in the 
following expression: 

a= lim l(E» &h -NE hulk ), (1) 

N—>oo Z 

where E^ ah is the total energy of an N-atom slab, i?buik 
is the total energy of the bulk per atom, and the fac- 
tor i accounts for the two surfaces in the slab unit-cell. 
However, this expression will diverge with increasing slab 
thickness if there are numerical differences between the 
calculation for the bulk and the slab (such as differences 
in the k-point mesh, etc.).— Two methods have been sug- 
gested to determine surface energies while avoiding this 
divergence. Boettger— uses the bulk energy in Eq. [T] 
as E^ h — E^^ x , thus avoiding a calculation on a sep- 
arate bulk system and effectively eliminating the errors 
from differences in k-point sampling, while Fiorentini and 
Methfessel^ 3 - make the assumption that in the limit of 
large TV one can rewrite Eq. [1] as: 

E^ h ^2<j + NE hulk . (2) 

If the total energy of the slab depends linearly on slab 
thickness N, the bulk energy term £"buik can be taken 
as the slope and used in Eq. [1] also avoiding a calcula- 
tion on a separate bulk system. In practice, it has been 
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TABLE II: Surface relaxations for the top three layers of all surfaces considered in this study. Reported are the values for 
a 13-layer slabs, compared to LDA and experimental values from the literature. For surfaces that are known to reconstruct, 
experimental values and those calculated as reconstructed in this study are given in square brackets. 
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shown that divergence can be avoided when large and 
matching k-point samplings are used for the slab and 
bulk calculations^ Here we compare these three meth- 
ods using the labels "Boettger" , "Fiorentini" , and "di- 
rect" respectively. Fig. [5]shows the comparison of surface 



energy versus slab thickness for Pd(100) surface energies. 
Here two different calculations are performed for the "di- 
rect" determination of -Ebuik- First, a single FCC unit 
cell with a 16 x 16 x 16 k-point mesh (dashed circles in 
Fig. ^ and second a two atom tetragonal unit cell with 
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FIG. 2: Surface energy versus slab thickness for the Pd(100) 
surface calculated with the methods of Boettges^, Fiorentini 
and Methfessel, and directly from Eq. [T]with two different 
values for i?buik- 
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FIG. 3: Surface energies of unrelaxed (dashed lines) and re- 
laxed (solid lines) slabs of Pd(110), Pd(100), and Pd(lll) 
surfaces versus slab thickness. 



a 16 x 16 x 11 k-point mesh to, as closely as possible, 
match that of the k-point mesh used for the slab (solid 
circles in Fig. [2|. 

It can be seen from Fig. [5] that the best non-divergent 
results as a function of slab thickness are obtained 
with the method of Fiorentini and Methfessel, while the 
method of Boettger oscillates around these converged val- 
ues. This periodicity is likely a manifestation of the finite 
size effects discussed earlier. Furthermore, even though a 
linear regression of the total energy versus slab thickness 
yields an excellent fit, the localized slopes (^iab — ^slab 1 ) 
of the Boettger method lead to energy fluctuations in the 
bulk reference energy of the order 1 x 10 -3 eV. These fluc- 
tuations are further amplified by the factor — N. These 
same fluctuations can be seen in all methods when com- 
paring surface energy data for any given slab thicknesses 
N and N ± 1; however, the choice of bulk reference in 
the direct and Fiorentini methods do not further amplify 



TABLE III: Surface energies for 13-layer slabs in both 
unrelaxed and fully relaxed geometries reported here in 
(eV/atom). LDA and experimental values compared with the 
surface energy of the relaxed surfaces reported here in units 
of (J/m 2 ). 
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them. Finally, neither the Boettger or Fiorentini meth- 
ods diverge in the same manner that the direct method 
does when the k-point grids do not match well for the 
bulk and surface calculations (dotted circles in Fig. [2]). 
However, similar to the findings of Da Silva et aL,— we 
find that if the same care is used when choosing the k- 
point mesh (as is the case for the solid circles in Fig. 
[2J an acceptable result is achievable through the direct 
method. For the remainder of this paper reported surface 
energies are calculated with the method of Fiorentini and 
Methfessel. 

Surface energies are calculated for both relaxed (but 
unreconstructed) and unrelaxed surfaces of all metals 
considered in this study. As an example, Fig. [3] shows 
the surface energies as a function of slab thickness for all 
three Pd surfaces, both unrelaxed (dashed lines) and re- 
laxed (solid lines). From the figure it is clear that surface 
energies are well converged for slabs as thin as 6-layers. 
Also, as expected, the surfaces that relax the most - i.e. 
(110) - experience the largest change in surface energy. 

Calculated values for 13-layer slabs for all metals stud- 
ied can be found in Table HTT1 All of the values calculated 
are converged with an accuracy comparable to those of 
Fig. [3J Furthermore, in the table we compare our present 
PBE calculated values with LDA and experimental val- 
ues from the literature. Remarkably, PBE seems to un- 
derestimate surface energies by as much as a factor of 2. 
If better comparison to experimental values is desired, 
LDA represents the better choice of exchange-correlation 
functional for all of the (111) surfaces for which data 
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was available. However, the converged value of surface 
energy for the Ti(0001) surface, calculated within PBE- 
GGA appears to more closely match the experimental 
data. Furthermore, regarding the (111) surfaces, our pse- 
dopotential calculations show excellent agreement with 
the FLAPW calculations- having only small differences 
all less than 0.03 eV/atom for the metals studied. 

In the context of poorly represented surface energies, 
we note the existence of exchange-correlation function- 
als that are designed to address this class of problem. 
First, one in which the surface and bulk of the material 
are treated with region specific functional with an in- 
terpolation region between them i 65 i 66 As Armiento and 
Mattsson point out, the functional they develop in their 
present article (AM05) give results similar to LDA, and 
while they have not been used on molecules, they rep- 
resent a good choice for solid state materials.— Second, 
in removing the bias toward the description of free-atom 
energies through a restoration of the second-order gra- 
dient expansion for exchange over a wide range of den- 
sities Perdew et al. create a functional built on PBE, 
PBEsol. 67 This same construction principle has more re- 
cently been applied to the meta-GGA functional (TPSS) 
by Perdew et al. and also shows promising results. 68 
The promise of these functionals (AM05 and PBEsol) 
regarding bulk properties is further supported by a re- 
cent overview of the accuracy of new functionals for bulk 
solids by Csonka et alf^ It is shown that both function- 
als perform notably better for most metals. Furthermore, 
regarding bulk Pd, which in our study shows the largest 
discrepancy with respect to experiment, they show an 
agreement of the lattice constant within 0.003 or less for 
both functionals. 69 Ropo et al. compare the two func- 
tionals for calculation of surface energies of metallic sur- 
faces, including the Pd(lll) surface, where they show 
an agreement with the experimental value of 2.00 J/m 2 
for both PBEsol (2.08 J/m 2 ) and AM05 (2.02 J/m 2 ). 70 
Further study of surfaces with these functionals are still 
needed, and the popularity of the functionals PBE and 
LDA still dominates the study of extended systems. 

D. Work Function 

As mentioned earlier, the work function is the mini- 
mum energy needed to remove an electron from the bulk 
of a material through a surface to a point outside the 
material, and can be written as: 

$ = K/acuum — Ep. (3) 

Again we need to test the calculation of the work function 
through use of the slab-supercell approximation. Thus, 
we compare and discuss the calculation of work func- 
tions a posteriori by two methods. First, we use directly 
Eq. [3l where both the potential in the vacuum region 
^vacuum and the Fermi energy Ep are derived from the 
same calculation. Second, we use a methodology utiliz- 
ing macroscopic averages ! 15 ' 71 Here we take the macro- 
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FIG. 4: A comparison of work function calculations methods 
for work functions of Pd(100) slabs versus slab thickness. The 
bulk reference used is a single unit cell of FCC Pd. 



scopic average of a potential step across the surface of 
the slab, AV = Kacuum - V s ub interior, and reference it to 
the macroscopic average of a potential obtained from a 
calculation on a separate bulk system, V^^ T °. This al- 
lows for the use of the Fermi energy from the bulk system 
^F,tmik, which will not suffer from quantum size effects.— 

* = AV + <fk r ° - £ F , bulk (4) 

Note here that the window for macroscopic averaging is 
the unrelaxed equilibrium layer spacing. All of the po- 
tentials discussed here (both from the slab and the bulk 
reference) refer to the electrostatic part of the total po- 
tential (the Hartree potential) . This part of the potential 
tails off more rapidly in the vacuum region of the slab- 
supercell when compared to the full Kohn-Sham poten- 
tial including the exchange-correlation potential. 

A comparison of these methods of calculation can be 
seen in Fig. [4] Here the work function of Pd(100) slabs 
is plotted versus slab thickness for both the unrelaxed 
slabs (dashed lines) and relaxed slabs (solid lines). It is 
clear from the figure that both methods converge quickly 
for the unrelaxed case (in as little as 7-layers). How- 
ever, when taking relaxation of all layers into account, the 
method of macroscopic averaging converges much more 
slowly. This is due to the fact that the averaging window 
for the reference potential and the layer spacing of the 
interior of the slab are not equivalent. Work functions re- 
ported in the remainder of this paper are taken directly 
by Eq. [3] without bulk references. 

As a demonstration of the convergence of the work 
function as a function of slab thickness, the results for 
three Pd surfaces are shown in Fig. Here it can be 
seen that the work function of the (111) and (100) sur- 
faces converges much more quickly than for the (110) 
surface. Even so, the work function of the (110) surface 
is converged to within O.leV by 8-layers and < 0.05eV 
by 11-layers. As in the previous subsection regarding the 
surface energy, oscillations can be seen. Furthermore, it 
is know that the work function of different facets can 
differ greatly^ 2 - and the general trend is that work func- 
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FIG. 5: Work function versus slab thickness for the (111), 
(100), and (110) Pd surfaces, both unrelaxed (dashed lines) 
and relaxed (solid lines). 



tion will decrease with decreasing layer packing, which is 
consistent with the findings in Fig. O 

Our final results for the work function of all the metals 
considered can be seen in Table HVl LDA and experimen- 
tal values from the literature are provided when available. 
With the exception of Au and Ti, LDA again represents 
the better choice of exchange-correlation functional if a 
better match to experimental results is desired. We find 
it particular interesting that Au and Pt should behave so 
differently with respect to the functional used. Regard- 
ing Al, we see the same trend as reported in Da Silva et 
al., where the work function of the (111) and (110) sur- 
faces are nearly equal, and the (100) has the larger of the 
three. Also, recall that the (100) and (110) surfaces of Au 
and Pt are likely in a reconstructed state for the experi- 
mental values. Finally, regarding the (111) surfaces, our 
psedopotential calculations again show excellent agree- 
ment with the FLAPW calculations^ having only small 
differences all less than 0.03 cV for the metals studied. 



IV. SUMMARY 



TABLE IV: Work functions calculated for 13-layer slabs of 
all metal faces considered in this study. The values in square 
braces are for surfaces that are likely reconstructed. All values 
in eV. 



Surface 


$ 


$LDA 


$Expt. 


Al (111) 


4.02 a ,4.04 d 


4.25 6 ,4.21 d 


4.23±0.02 c 


(100) 


4.30 a 


4.38 6 


4.42±0.03 c 


(110) 


4.09° 


4.30'' 


4.12±0.02 c 


Pd (111) 


5 25™ 5 22 d 


5.64 d 


5 90+0 01 m 5 55+0 01 c 


(100) 


5.11° 




5.65+0.01™ 


(110) 


4.87° 




5.20+0.01™ 


Pt (111) 


5.69 ,5.69 d 


6.06 d 


6.08±0.15*,6.10±0.06 / 


(100) 


5.66° 




[5.82+0. 15 fc ] 


(110) 


5.26 a 


5.52 9 


[5.35+0.05'] 


Au (111) 


5.15" 


5.63 h 


5.26+0.04* 


(100) 


5.10" 


5.53 ft 


[5.22+0.04*] 


(110) 


5.04" 


5.41 fe , 5.39 9 


[5.20+0.04*] 


Ti (0001) 


4.38 a ,4.40 d 


4.66 d 


4.33 J 


"Present study 



6 PWPP-LDA-SlabS 
c Ref^ 

d FLAPW-LDA-Slab 10 

e Ref-£ 

/Refl£ 

»PWPP-LDA-Slab£i 
h PWPP-LDA-Slab2£ 
^Refl 8 . 

JRcf 79 
fcRcf 80 
( Rcf 81 
m Rcf _82 



well quantitatively to experimental values. This may be 
understandable for the surface relaxations and surface 
energies, where even the experimental values have some 
room for interpretation. While for the work functions, a 
property for which the experimental values are more re- 
liable, we have found that neither LDA or GGA is unan- 
imously a better choice. 



Calculations of the relaxations, surface energies, and 
work functions of low index metallic surfaces were made 
for the (100), (110), and (111) surfaces of Al, Pd, Pt, 
and Au and the (0001) surface of Ti using finite slab 
approximations. We have paid particular attention to 
the issues that arise from slab thickness and consequently 
finite size effects, and we have shown that convergence 
of these quantities can be achieved for slab thicknesses 
that are greater than 6-layers for (111), (100), and (0001) 
surfaces and 10-layers for (110) surfaces. We find that the 
use of bulk references for calculations of surface energies 
and work functions can be detrimental to convergence, 
especially when surface relaxations are being considered. 

We compared our results within the GGA with exper- 
imental and LDA results from the literature. Our results 
and comparison have shown that, even though converged 
values can be achieved, calculated values do not match 
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APPENDIX: MAGNETIC BULK PALLADIUM 

Recent calculations in the literature have shown that 
when evaluated within the generalized gradient approx- 
imation (GGA) Pd yields an erroneous magnetic bulk 
ground state i 25 ' 26 In light of these observations we have 
performed a similar evaluation of the ground state of Pd 
using PBE-GGA and the pseudopotential mentioned in 
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the text of this paper. In a similar manner to Alexan- 
dre et al. we find also a magnetic ground state for bulk 
Pd. We plot both the total energy and the total mag- 
netization as a function of FCC lattice parameter for Pd 
evaluated with GGA in Fig [6] and for comparison the 
same plot is generated for Pd within LDA in Fig [JJ 



Lattice Constant (A) 



FIG. 6: Total energy per unit cell and total magnetization 
per unit cell plotted versus the FCC lattice constants of bulk 
Pd using the GGA-PBE exchange-correlation functional with 
spin-polarization. 




Lattice Constant (A) 



FIG. 7: Total energy per unit cell and total magnetization 
per unit cell plotted versus the FCC lattice constants of bulk 
Pd using the LDA exchange-correlation functional with spin- 
polarization. 



Finally, we have calculated the surface properties of 
the Pd(100) with spin-polarization. In all cases the slab 
maintains a total magnetization, however, this has little 
effect on the converged surface properties discussed in 
the body of this paper. The changes in the relaxations 
and surfaces energies are negligible, while the change in 
the work functions is seen as an increase of 0.5%. 
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